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Exploring two non-perturbative definitions of ca*^ 

S. Dmr a , M. Delia Morte a 

a DESY, Platanenallee 6, 15738 Zeuthen, Germany 

We present two determinations of the coefficient ca in quenched QCD, needed to build the O(a) improved axial 
current. The first condition used is the requirement that the PCAC quark mass, as a function of Xq, stays flat for 
a non-trivial spatial phase for the fermions in the Schrodinger functional. The second condition is that the PCAC 
relation for the ground-state and the first excited state at finite L give the same quark mass. Our results confirm 
previous findings that in the quenched theory the intrinsic 0(a) ambiguity of ca gets relevant around /3 ~ 6.0. 



1. INTRODUCTION 



The improved axial current 

4(x 



A (x)+ac A ^^P(x) 



(1) 



with Aq(x) =-0^7075^ and P(x) =^75^ is de- 
signed to modify the scaling behavior of on-shell 
quantities from 0{a) to 0(a 2 ), if ca and csw 
are chosen appropriately. Therefore, an accurate 
determination of the coefficient ca is an impor- 
tant ingredient in the improvement programme 
a la Symanzik. Here, we present results in the 
quenched theory, but the real motivation is to 
find a criterion which is practical in a dynami- 
cal setting, where large cut-off effects have been 
found [J. 

Our data have been generated in the Schro- 
dinger functional (SF) setup: we use a TxL 3 box, 
applying some Dirichlet type boundary condi- 
tions in the time direction (i.e. at xq — 0, T) while 
keeping the gauge-field periodic in space and the 
fermions periodic up to a phase: U^Xq, x+Le^) = 
U^x), ip^(x ,x+Le^) = exp(i6>)Vv(a;). The ini- 
tial pion with flavor content ij is created at xq — 
through the boundary operator 



a 6 

°« = 7s z]Ci(u)75Ci(v) w(u-v) 



L 3 

U,V 



(2) 



introducing the wavefunction u> w.r.t. the rela- 
tive position of the boundary fields Q , Q , and ab- 
sorbed by the local current Aq(x) in the bulk or 
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by the analogous operator Oy at xq — T. This 
means that we consider the SF correlators 



fx(x ,T,L) = -L.(X(x)0) 



h(T,L) 



- 2 (o'o) 



(3) 
(4) 



with X cither Aq or P to get the PCAC mass 
m = r- 



r(x ) 
s(x ) 



- acA s with 

2f P (x ) 
d Q d^fp(x Q ) 



(5) 
(0) 



2f P (x ) 

Requiring m constant at fixed /3 leads to a def- 
inition of the improvement coefficient through 
c,4 = Ar/As, where the difference may be w.r.t. 

(0) 6 at fixed x =T/2 (ALPHA [5] ) 

(1) xo at fixed 9 (LANL [3], "slope criterion") 

(2) state in 0[w] (UKQCD 4, "gap criterion" ) 

and depending on the choice, physical quantities 
differ by 0(a 2 ) effects. This means that there 
is an intrinsic 0{a) ambiguity in ca itself which, 
already for /3~6, is not such a small effect [314] . 

Here we investigate (in a quenched setting) 
which one, out of (0)-(2), might be a promising 
criterion for N{ = 2, with a view on the following 
wishlist: (i) no high-energy state involved, i.e. no 

< r o (say) used, (ii) large "sensitivity" , i.e. not 
too small value of As, (iii) affordable numerical 
effort, i.e. not requiring large volume, (iv) "scal- 
ability" , i.e. allowing to move to another /3 while 
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Figure 1. r(cco) and s(xrj) from a SF wall source with four f?-values at /? = 6.2, re = 0. 13485. 



keeping physics in units of ro constant. Both the 
ambiguity of improvement coefficients and how to 
deal with it have been discussed in [5]. 

We emphasize that the SF states generated by 
the boundary operator @ are multiplicatively 
renormalizable. The associate Z-factor cancels 
in the ratios JSHHJ), and everything is scalable. 

2. THETA CRITERION 

The old ALPHA criterion [2] resulted in rather 
small As values. Furthermore, for Nf — 2 several 
6>-angles mean several simulations. Therefore, we 
didn't investigate (0) further. 

3. SLOPE CRITERION 

The slope criterion (1) requires only one re- 
value, but for completeness we decided to test it 
for (9 = 0,7r/3,27r/3,7r. 

Fig. 1 shows r(xo) and s(xq) in the SF. For 

(3=6.2, 12 3 x27, wall source 
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Figure 2. -ca via slope in Fig. 1 with x re f = 15. 



9 = ir there is a good sensitivity. Using Ar(xo) = 
r(xo)-^r(x vc {) with x rc f around the extremum (and 
ditto for As) the recipe ca = Ar/As yields a lo- 
cal ca{xo). The remnant dependence on x rc f was 
checked to be small. This and the dependence on 
0, L,T, k,u> represent genuine 0(a) effects on ca- 
Fig. 2 displays —ca(xq) determined via this 
"slope criterion" . For 9 = n and small enough 
L there is an early plateau which is a sign that 
all states but the lowest two have disappeared. 
The corresponding value for ca seems consistent 
with the old ALPHA determination [2]. 

4. GAP CRITERION 

From a transfer matrix analysis one gets 

f x ~^pZe- M ^{l+r,le- Ax °+r? x e- M °^} 

f^\ P 2 e- M * T (7) 

with the matrix element £ = (0, 0| J&T|7r, 0) and 
known representations of p, rj x ° in terms of states 
\Q,n) with a given set of quantum numbers and 
excitation level. In A denotes the gap in the 
corresponding (here: pseudoscalar) channel and 
Mq the mass of the lowest (0 ++ ) glucball state. 
The important point is that the coefficients p, rj x 
depend on the initial state and hence on the 
wave function u>, while £,?7x °-0 not. This cre- 
ates the possibility to linearly combine the cor- 
relators (JJJ over several u> to build one which is 
clearly dominated by either the ground-state or 
the first excited state. After checking that the 
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(3=6.0, 12 x12, 4 WFs and projection to grd or exc state 
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Figure 3. r(xo) and s(xq) from the ground-state and the first excited state at /3 = 6.0, k = 0. 13415, 9 = 0. 



corresponding effective masses are indeed distinct 
(i.e. Ma < M*), one may define Ar as the differ- 
ence of the expressions 10 w.r.t. these two linear 
combinations (at a given xq, 6) and ditto for As. 
This amounts to an operational definition - time 
slice by time-slice - of ca according to the "gap 
criterion" (2). 

Fig. 3 shows r and s at 9 = for the ground- 
and excited state, built (a posteriori) from the 
four hydrogen-type wave functions used in the 
simulation. It seems this method yields a rather 
good sensitivity As. In the supposed ground- 
state both r(xo) and s(xq) happen to be almost 
flat from a rather early time on. 

Fig. 4 shows — ca(xq) determined via this "gap 
criterion" . The plateau is reached at xo ~ 0.75 fm, 
and typical values are somewhat smaller than the 
original ALPHA value at f3 = 6.0 0. The absolute 
value \ca\ — —ca that we get decreases with L, 
but for L > 3ro this effect seems not particularly 
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pronounced any more. We take this as a sign that 
there is good hope to determine ca in a (1.5 fm) 4 
box in the N{ = 2 theory. 

5. SUMMARY 

We have tested several improvement condi- 
tions that might be used to determine ca non- 
perturbatively in an unquenched setting. The 
0(a) ambiguity of ca that has been pointed out 
previously |3I4| has been confirmed via, (a) ex- 
tension of (1) to several 9 and, (b) refinement 
of (2) with a wave function projection technique 
which achieves dominance by ground- or 1st ex- 
cited state in the region 0.5 fm ... 1.0 fm. We hope 
this proves sufficient to determine ca with Nf = 2. 
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Figure 4. -ca via gap in Fig. 3; ditto at L — 16, 24. 



